Shot noise in tunneling transport through molecules and quantum dots 
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We consider electrical transport through single molecules coupled to metal electrodes via tunneling 
barriers. Approximating the molecule by the Anderson impurity model as the simplest model which 
includes Coulomb charging effects, we extend the "orthodox" theory to expand current and shot 
noise systematically order by order in the tunnel couplings. In particular, we show that a combined 
measurement of current and shot noise reveals detailed information of the system even in the weak- 
coupling limit, such as the ratio of the tunnel-coupling strengths of the molecule to the left and 
right electrode, and the presence of the Coulomb charging energy. Our analysis holds for single-level 
quantum dots as well. 

PACS numbers: 73.63.-b, 73.23.Hk, 72.70.-|-m 



I. INTRODUCTION 

Single-molecule electronics provides an exciting pos- 
sibility for further miniaturization of electronic devices. 
Nonlinear electrical transport through single molecules 
deposited between mechanically controlled break junc- 
tions has been measured in several experiments4iS*i^ 
Furthermore, recently observed diode behaviori^ mem- 
ory effects and negative differential conductance^ in de- 
vices which make use of molecular films have attracted 
much interest. 

A crucial problem for understanding and designing 
molecular electronic devices is the lack of control of 
the metal-molecule interfaces. In molecular-film exper- 
iments, deposition of the covering electrode is very diffi- 
cult and can easily destroy the film or induce short cuts. 
Once the molecules are deposited, no direct study and 
manipulation of the contact is possible. Even in a break 
junction experiment one has very limited influence on 
how a molecule bonds to the two electrodes. Information 
on the coupling strengths rL,rR between the molecule 
and electrodes has to rely on transport measurements. 
From the measurement of the current /, however, only a 
combination of Fl and Fr can be deduced, whereas the 
asymmetry ratio Fl/Fr remains undetermined. 

In addition, it is impossible to judge from the I ~ V 
curve alone, whether charging effects such as Coulomb 
blockade are relevant to the transport. This is because 
any observed set of current steps expected for an / — 
curve displaying Coulomb blockade can be mimicked by 
transport through a set of noninteracting levels. As so 
far there is no quantitative theory for describing hybrid 
metal-molecule functions, it remains a matter of opinion, 
whether such a fit is reasonable for a given experiment. 

Further information can be gained from the shot noise 
S, allowing, in principle, to determine the ratio of 
the couplings. Furthermore, the combination of cur- 
rent and shot-noise measurement provides information 
about the presence and size of Coulomb charging en- 
ergy. For transport through a noninteracting levelSi2iili 
the so-called Fano factor F = S/2eI is known to be 



(F^ + r|)/(FL + Fr)2. However, transport through 
weakly coupled molecules can not be modeled by non- 
interacting molecular levels, since Coulomb blockade ef- 
fects become important 

Some results have been obtained before for shot noise 
through interacting systemsii^ii^*iSiiLi& However, the fo- 
cus in these works was on transport through impurity 
states or quantum dots, where some questions pertinent 
to molecules are irrelevant. 

In the present paper, we reformulate the "orthodox" 
theory of shot noiseiS in a way that allows a system- 
atic perturbation expansion in the coupling strength. We 
study in detail the shot noise for a molecule accommo- 
dating one level, described by the single-level Anderson 
model, in the lowest-order perturbation theory. We per- 
form a complete classification of all transport regimes, 
and evaluate the corresponding Fano factors, recovering 
known expressions for the cases previously studied in the 
literature. The Fano factor turns out to depend in a 
nonmonotonic way on the asymmetry ratio Fl/Fr. The 
Fl/Fr dependence of the Fano factor can be used to 
identify the different transport regimes. Our results hold 
for transport through single-level quantum dots as well. 
In contrast to molecular devices, the contacts to leads 
in the semiconductor quantum dots are usually well con- 
trollable. 



II. THE MODEL 



The simplest model for electron transport through 
a molecule with Coulomb interaction is the Anderson- 
impurity model described by the Hamiltonian H = Hj^ + 
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Hn + Hmo\ + Ht.l + Ht.b. with 

Hr = ^'^ekcrrdl^r'^kiTr, (1) 
kcr 

Hmo\ = y^£cr4cff + Un-\ni, (2) 

(J 

HT,r ='^{tral^^Ca + h.C.) (3) 
ka 

and r = L, R. Here, Hl and ifft niodel the left and right 
electrodes with noninteracting electrons, Hmo\ describes 
the molecule with one relevant (spin-dependent) molec- 
ular level of energy e^, and the Coulomb interaction U 
on the molecule {n-\ and are the number operators 
for electrons with corresponding spin). Tunneling be- 
tween leads and molecule level is modeled by i?T,L and 
Ht.b.- The coupling strength is characterized by the in- 
trinsic linewidth = 27rpe|ij.p, where Pe is the (con- 
stant) density of states of the leads. Furthermore, we 
define the total linewidth F = Fl -I- Fr. The Fermi oper- 
ators flfco-rC'^fef'") ^^'^ cJ.(C(t) create (annihilate) electrons 
in the electrodes and the molecule. 

We are interested in transport through the molecule, in 
particular in the current / and the (zero-frequency) cur- 
rent noise S. They are related to the current operator I = 
(-^r-^l)/2, with Ir = -i{e/h) J^ko- (^'■"Lr'^'^ ~ ^■'^■) hie- 
ing the current operator for electrons tunneling into lead 
r, by / = (/), and 

/oo 
dt{si{t)di{o) + si{o)di{t)) , (4) 
-oo 

where Si{t) = i{t) - (/). 

III. DIAGRAMMATIC TECHNIQUE 

A diagrammatic technique to perform a systematic 
perturbation expansion of the current through localized 
levels has been developed in Ref. Here, we expand 
this technique to address the perturbation expansion of 
the current noise as well. After this general formulation, 
we specify the results for lowest order in F and discuss im- 
plications for transport through molecules in this regime. 

All transport properties are governed by the nonequi- 
librium time evolution of the density matrix. The leads 
are described as reservoirs of noninteracting electrons, 
which remain in thermal equilibrium. A finite transport 
voltage enters the difference of the electrochemical po- 
tentials /ir for the left and right leads, r = L, R. These 
degrees of freedom can be integrated out making use of 
Wick's theorem, such that the Fermion operators are con- 
tracted in pairs. As a result, we end up with a reduced 
density matrix for the molecule's degrees of freedom only, 
labeled by x- In our case, four molecule states are pos- 
sible: the level is empty, x — singly occupied with 
either spin up or down, x —1 or |, or doubly occupied 



X = d. The time evolution of the reduced density matrix, 
described by the propagator H^/^(i',t) for the propaga- 
tion from state x at time t to state x' at time i', can be 
visualized by diagrams. 
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FIG. 1: An example for the time evolution of the reduced den- 
sity matrix. The upper and lower lines represent the forward 
and backward time propagation along the Keldysh contour, 
respectively. Tunneling lines correspond to the reservoirs L, R 
connecting pairs of vertices. The resulting changes between 
the four molecular states are indicated. 

An example is shown in Fig. ^ Vertices represent- 
ing tunneling, Ht, are connected in pairs by tunneling 
lines. The full propagation is expressed as a sequence of 
irreducible blocks W^i^(t' , t) containing one or more tun- 
neling lines. They are associated with transitions from 
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FIG. 2: An example for an irreducible block W. 

state X at time t to state %' at time t' (see Fig. 13) • This 
leads to the Dyson equation for the propagator 

n(0,i) = l+ / dt2 dtiW{t2,ti)U{ti,t), (5) 

where the boldface indicates matrix notation related to 
the molecular state labels. To calculate transport prop- 
erties in the stationary limit, we start at time t — > — oo, 
at which the system has a given diagonal but other- 
wise arbitrary distribution of probabilities p^'' to be in 
state Xi comprised in the vector p'"'*. The stationary 
probability distribution p*^' = 11(0, —oo)p'"'*, however, 
does not depend on p™'*, i.e., H^'^(0, — oo) = p^,, in- 
dependent of X- To determine p''' we use the Dyson 
equation Eq. ^ in the limit t — oo, multiply p'"'* 
from the right, and define W — J_^dtKW{0,t) to get 
pst _ pinit _ iWpSt^ where 77 ^ 0+ comes from a 

convergence factor exp(— ?7|f2|/^) in the time integral. 
This leads to Wp*^* = together with the normaliza- 
tion e-^p*^* — 1, where the vector e is given by = 1 
for all X- The matrix W has zero determinant and hence 
cannot be inverted (this can be seen from the sum rule 
e^W = 0). To write down a single matrix equation 
which determines p^* we make use of the normalization 
of probabilities and introduce the matrix W, which is 
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identical to W but with one (arbitrarily chosen) row xo 
being replaced with (F, ...,r). The stationary probabili- 
ties are the solution of the linear equation 



(Wp^ 



(6) 



This can be solved for p*** by inverting W. 

For a well-defined perturbation expansion in powers k 
of the coupling strength F we write W ~ ^kLi'^^'^K 
W = W('=', and p"' = X;^o P'*'''''- The order in 

F corresponds to the number of tunnel lines contained 
in the irreducible blocks. As a consequence, W and W 
start with first order in F, whereas p''* starts in zeroth 
order. The zeroth-order stationary probabilities are 



„st(0) 



F 

'XX(1 



(7) 



and higher-order corrections are obtained iteratively by 



st(fc) 



m=0 



m+l)pSt(m) 



(8) 



for k — 1,2,.... The diagrammatic rules for the calcula- 
tion of the blocks W'^'^' will be summarized at the end of 
this section. 

The formulation of the perturbation expansion of the 
stationary probabilities is the first important step for a 
perturbation expansion of the current. The diagram- 
matic representation of contributions to the current in- 
volve a block W^, in which one (internal) vertex is re- 
placed by an external vertex for the current operator / 
divided by e/h. The current can also be expanded order 
by order in the coupling strength F, 



fc-i 



2h 



m=0 



(9) 



for k = 1,2, . . ., where the factor 1/2 corrects for double 
counting of the external vertex being on the upper and 
lower branch of the Keldysh contour. 

Now, we turn to the evaluation of the current noise 
5 = 2/°^ dt [{I{t)i{0)+I{0)I{t))-2{I)% which involves 
expectation values of two current operators at different 
times. These two current operators can either appear in 
the same irreducible block, denoted by W^^, or in two 
different blocks. In total, the current noise is expressed 
as 

S = j-e^ (W^^p^* + W^PW^p"*) , (10) 

(where the factor 2 in the above definition has canceled 
against the correction factor 1/2 for double counting) 



with 



rfti [n(o,<) -n(o,-oo)] . (11) 

The constant part 11(0, — oo) in the integrand is associ- 
ated with the (7)^ contribution to the noise. Further- 
more, it guarantees that the integral is convergent and 
that P exists. 

We also can formulate a perturbation expansion of P 
in orders of F. To do so, we first relate P to the blocks 
W and the stationary probabilities p^* by 



WP 



(12) 



which follows from the Dyson equation and Wp'' = 0. 
As discussed above, W is not invertible, i.e., in addition 
to Eq. I|12|) we need some extra condition to determine P. 
This extra condition is e-^P = 0, which follows from the 
definition, Eq. (|11() . together with the Dyson equation, 
Eq. ||SJ|, and e-^W = 0. Introducing the matrix Q with 
elements given by 



Qx'x = bx' -'^x',x)(l -<5x',xo), 



(13) 



Eq. (|12|l can be recast as WP = Q, which allows the 
determination of P by inversion of the matrix W intro- 
duced above. The solution for P can be expanded in 
orders of F. We observe that, since the expansion of Q 
starts in zeroth order in F, P starts in order F^^, 



>(-i) 



(W(i')-iQ 



(0) 



(14) 



This reflects the fact that between the external cur- 
rent vertices the system evolves in time, as described 
by n(0,f). The time scale at which the dot relaxes 
and correlations decay is set by the inverse of the tunnel 
rates, i.e., by h/T. As a consequence, the time integral in 
Eq. (|ll|l is effectively cut at this time scale, i.e., the dom- 
inant contribution to P is of order F~^. Furthermore, it 
is crucial that P starts in order F^^, as otherwise the 
second term in Eq. (|10|) would not contribute to noise to 
first order in F (cf. Ref. 20), leading to a wrong result 
even for non- interacting systems {U — 0), where exact 
expressions are available for comparison, 

Corrections of higher order, k — 0, 1, can be com- 
puted from 



Q(fc+1) _ ^ ^^[k-m+l)-p 



(m) 



(15) 

Eventually, we have all ingredients at hand to expand the 
noise in orders of F, and we find 
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k — 1 I k — rnk — m — m' — l 

m' — l m" — — 1 



(16) 



r 



for A; = 1, 2, . . .. In lowest order this result simplifies to 
5(1) = (eVS)e^ (W"W +W^(i)p(-i)W^(i))p^*(0). 

In summary of the technical part, we have derived a 
scheme to perform a systematic perturbation expansion 
of the current and the current noise. First, one calculates 
the irreducible blocks W, W^, and W^^. Second, the 
stationary probabilities p*^' are obtained from Eqs. {Tj) 
and ^ and P follows from Eqs. ifTHl and lfTK|l . Eventu- 
ally, current and current noise are given by Eqs. © and 

The irreducible blocks W are calculated according to 
the rules specified in Rcf. 19 [where the irreducible blocks 
were denoted by S and included a factor i in the defi- 
nition, S = iW]. We complete this section with a sum- 
mary of these diagrammatic rules, given as follows: 

(1) For a given order k draw all topologically different 
diagrams with 2k vertices connected by k tunneling lines 
(an example of a first- and second-order diagram is shown 
in Fig.jSJ. Assign the energies to the propagators and 
energies lui {I = 1, k) to the tunneling lines. 

(2) For each of the {2k — 1) segments enclosed by two 
adjacent vertices there is a resolvent l/{AEj+iO^), with 
j — 1, ...,2k — 1, where AEj is the difference of the left- 
going minus the right-going energies. 

(3) The contribution of a tunneling line of reservoir r is 
{l/2'K)Trf{uJi — /ir) ff the line is going backward with re- 
spect to the closed time path and (l/27r)r,.[l — /(w/ — /i^)] 
if it is going forward. Here, f{x) is the Fermi function. 

(4) There is an overall prefactor (— i)(— 1)^, where c is 
the total number of vertices on the backward propagator 
plus the number of crossings of tunneling lines plus the 
number of vertices connecting the state d with "f . 

(5) Integrate over the energies ui of the tunneling lines 
and sum over all reservoir and spin indices. 

The blocks and W-^^ are determined in a similar 
way. The only difference to W is that in W'^ (W^^) one 
(two) internal vertices are replaced by external ones rep- 
resenting Ih/e. This amounts to multiplying an overall 
prefactor, which arises due to the definition of the cur- 
rent operator and since the number of internal vertices 
on the backward propagator may have changed. We get a 
factor +1/2 for each external vertex on the upper (lower) 
branch of the Keldysh contour which describes tunneling 
of an electron into the right (left) or out of the left (right) 
lead, and —1/2 in the other four cases. Finally, we have 
to sum up all the factors for each possibility to replace 
one (two) internal vertices by external ones. 
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FIG. 3: Current I and shot noise 5 vs. voltage for T — lOOK, 
ei = O.leV, = 0.2eV, U = 0.4eV and Tl/Fr = 1/1 (solid 
line), 1/3 (dashed line), and 3/1 (dotted line). The height 
of the plateaus labeled by z = 0, ..,4 are given in Table 
The curves are normalized to /max = {e/h)T /2 and Smax = 
(e^/fi)r/2, respectively. 



IV. RESULTS 

In the following we discuss shot noise in the lowest- 
(first-) order perturbation theory in F. We consider 
transport through a single level and allow for a finite spin 
splitting. Such a spin splitting could be realized due to 
the Zeeman effect of an external magnetic field (for semi- 
conductor quantum dots) or an intrinsic exchange field 
due to atoms with magnetic moments (for molecules). 
The molecule can acquire four possible states: the level 
being unoccupied, occupied with either spin t or |, or 
doubly occupied. The molecular level is characterized 
by level energies and e|, the Coulomb repulsion J7, 
and the coupling strengths Fl and Fr to the electrodes. 
In addition, the electron distributions in the metal elec- 
trodes are governed by electrochemical potentials /il and 
/iR, and temperature T . We choose a symmetric bias volt- 
age, such that = eV/2 and /xr = —eV/2. In this 
simple model without real spatial extent, the voltage is 
dropped entirely at the electrode-molecule tunnel junc- 
tions, meaning that the energies of the molecular states 
are independent of the applied voltage even if the cou- 
plings are not symmetric. For asymmetric coupling this 
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might not be entirely realistic. However, we do not wish 
to include a heuristic parameter to describe the possibly 
unequal drop of the bias. In a comparison to an actual 
experiment this could be easily amended. We also do not 
include an explicit gate voltage, as its effects are straight- 
forward to anticipate. 

In Fig. Owe plot the current and the current noise for a 
special choice of system parameters as a function of trans- 
port voltage. For this figure we consider a set of energy 
parameters (e^ = 0.1eV,e| = 0.2eV, U = 0.4eV) that is 
comparable to energies encountered in small molecules. 
A rescaling of energies by a factor 1/100 would be neces- 
sary to obtain energies of the order achievable in semicon- 
ductor quantum dots and by external magnetic fields. If 
the same scaling is applied to the couplings and the tem- 
perature, the figure would remain unchanged. This is a 
special feature of first-order transport. 

Electron transport becomes possible when charge ex- 
citations on the molecule become energetically allowed. 
Generally, at low bias, transport is exponentially sup- 
pressed, unless a degeneracy of states with different net 
charge is present (this could be tailored by application of 
a gate voltage, which we do not consider here). Each time 
when one of the four excitation energies ej^, e|, -I- J7, or 

+ U enters the energy window defined by the electro- 



chemical potentials of the electrodes, a transport channel 
opens. This gives rise to plateaus, separated by thermally 
broadened steps. Due to the symmetric application of the 
bias, the steps occur at voltages of twice the correspond- 
ing excitation energy. The plateau heights depend on 
the coupling parameters Fl and Fr only, i.e., they are 
independent of U and T. The analytic expressions for 
current, noise, and Fano factor of these plateaus, which 
are labeled by i = 0, 4, are given in Table |l| 

The curves in Fig.Oare normalized to /max = (6/^)^/2 
and iSmax = (e^/fi)r/2, respectively, which is reached in 
the large-bias limit for symmetric coupling Fl = Fr . For 
asymmetric coupling, the plateaus are reduced in height. 
In Fig. 13 wc show the results for 3Fl = Fr (dashed lines) 
and Fl = 3Fr (dotted lines) together with the case of 
symmetric coupling. The symmetry of our setup implies 
that all plateau heights are invariant under simultane- 
ous exchange of Fl with Fr and /j,l with /iR. However, 
the plateau height can change if only Fl and Fr are ex- 
changed, or if only the bias voltage is reversed, as shown 
in our example for the two plateaus labeled by 2 and 
3. This opens the possibility to access the asymmetry 
Fl/Fr experimentally by reversing the bias voltage and 
comparing the plateau heights. 
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TABLE I: Current, shot noise and Fano factor for the different plateaus in the current-voltage characteristic shown in Fig. El 
The plateau values depend only on the coupling parameters Fl.r. 



Figure 13 shows the result for one special choice of energy 
parameters and corresponding excitation energies. Nev- 
ertheless, Table|2is complete in the sense that it contains 
all possible plateau values for any configuration of the ex- 
citation energies relative to the electrochemical potentials 
of the electrodes. The classification of the configurations 
and the algorithm to find the corresponding analytic ex- 
pressions in Table is given in Table ^ Without loss 
of generality we restrict ourselves to J7 > and e-\ > . 
The different configurations are classified by specifying 
which excitation energies lie within the energy window 
defined by the chemical potentials /LtL,A*R. We find 13 



different possibilities, as listed in Table Ull For each case, 
the index i indicates the column where the correspond- 
ing analytic expressions for current, noise and Fano factor 
can be found in Table The indices 2* and 3* refer to 
columns 2 and 3 but with Fl and Fr being exchanged. 

In order to illustrate the use of the table we sketch 
the situation 2a in Fig. 0] realized in Fig. |31 in the re- 
gion between 0.4V and IV. In this situation transport 
through both spin states is present as the molecule gets 
charged/uncharged in the sequential tunneling events. 
Double occupancy is still out of reach, since the exci- 
tation energies -t- 1/ and + U are outside the energy 
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TABLE 11: Classification of all possible configurations that 
are possible for U > and et > ej,. A cross (x) or minus ( — ) 
indicates that the corresponding excitation energy lies within 
or outside the energy window defined by the electrochemical 
potentials /xl and fin, respectively. For each configuration, 
the analytic expression for current, noise and Fano factor can 
be found in Table in column i. The indices 2* and 3* refer 
to column 2 and 3 with Fl and Fr being exchanged. 



window opened by the applied voltage. 

8t+ U 

El+U 

^L= eV/2 

et 

£i 



--- I^R=-eV/2 

FIG. 4: A sketch of the configuration 2a listed in Table 
The excitation energies and e| lie in the energy window 
defined by the electrochemical potentials fj,h and /xr, and the 
energies + U and e^- + U lie outside. 

We now turn to the discussion of the Fano factor 
F = S/2eI. In Fig. [S] we show the Fano factor as a 
function of bias voltage for the same parameters as in 
Fig. 13 Again we show three curves with Fl /Fr = 1/1 
(solid line), 1/3 (dashed line) and 3/1 (dotted line). At 
small bias, eV <ti ksT, the noise is dominated by ther- 



mal noise, described by the well-known hyperbolic cotan- 
gent behavior which leads to a divergence of the Fano 
factorii2ii& The plateau for bias voltages below 0.2V cor- 
responds to the Coulomb-blockade regime, where trans- 
port is exponentially suppressed (case in Table ITU) . In 
the region between 0.2V and 0.4V (case la) transport 
through only one spin state (|) is possible. The Fano 
factor Fl for this case has been derived earlier in Ref. 0. 
For very large bias (region F4), all states of the molecu- 
lar level are involved in transport, and the Fano factor F4 
is again identical to the well-knownifl formula for trans- 
port through a resonant level in the absence of Coulomb 
charging energy. 
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FIG. 5: Fano factor vs. bias voltage for for the same param- 
eters as in Fig. 01 namely T = lOOK, = Q.leV,e| = 0.2eV, 
U = 0.4eV, and Fl/Fr = 1/1 (solid line), 1/3 (dashed line), 
and 3/1 (dotted line). The labels 0, la, 2a, 3a, 4 refer to the 
cases listed in Table |IT| 

For the regions in between, corresponding to the cases 
2a and 3a, the Fano factor is different. The expression 
for F2 has been very recently derived in Ref. llTL while F3, 
corresponding to region 3a, is, to the best of knowledge, 
a new result. 

On similar grounds as for the current and the noise, we 
find that the expressions for the Fano factors F2 and F3 
are not invariant under exchange of Fl and Fr alone, and 
also not invariant under reverse of the bias voltage alone. 
This is clearly seen in Fig. [S] in the very different plateau 
heights of the dotted and dashed curves, corresponding 
to exchange of Fl and Fr. Furthermore, we see that all 
plateau heights lie between 1/2 and 1, and that the Fano 
factor, in general, is a nonmonotonic function of the bias 
voltage. We find that Fi = F4 and F3 > F2 always hold, 
whereas Fi > F2 for Fl/Fr < 1/\/5 and Fi > F3 for 
Tl/Fr < 1/2 only. 

As a consequence, the pattern of the plateau sequence, 
in particular the relative height of the plateaus F2 and 
F3 compared to Fi = F4, indicates not only the presence 
of an interaction or charging energy, but also gives over- 
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complete information on the ratio Fl /Fr of the couphng 
strengths. This could be used in experiments to deter- 
mine these parameters in a consistent way. In an experi- 
ment observing more than one plateau, the overcomplete- 
ness would give narrow constraints (due to experimental 
uncertainty) on whether a single interacting level can ex- 
plain the measured values. Of course, it is always possible 
to fit n plateaus with n noninteracting levels and differ- 
ent couplings per level, so an absolute decision on the 
presence of interactions is not possible without the addi- 
tional application of a gate voltage. However, if a fit with 
an interacting level is feasible, the principle of parsimony 
should favor the model with fewer parameters. 




bias voltage [V] 



FIG. 6: Fano factor vs. bias voltage in the absence of spin 
splitting. We chose T = lOOK, = = O.leV, U = 0.4eV, 
and Fl/Fr = 1/1 (solid line), 1/3 (dashed line), and 3/1 (dot- 
ted line) . Also shown are the two cases for the corresponding 
noninteracting system {U ~ 0). 

In the absence of a spin splitting, e-f = e^, the number 
of plateaus is reduced. An example for the Fano factor 
as a function of bias voltage is displayed in Fig. |H1 The 
solid curve is for symmetric, and the dotted and dashed 



curves are for asymmetric coupling. Again, a nonmono- 
tonic sequence of plateau heights, F4 > F2, indicates the 
presence of Coulomb charging and Fl/Fr < 1/V2. For 
comparison we remark that for both negligible spin split- 
ting and zero charging energy (U = 0), only the plateaus 
_Fo and F4 appear, the plateau heights are invariant under 
reverse of bias, and no nonmonotonic behavior is seen, as 
also shown in Fig. El 

Finally, we briefly comment on the steps connecting 
the plateaus of the Fano factor. In the present first- 
order approximation the steps are solely broadened by 
temperature. In addition, we observe that sometimes a 
peak in the Fano factor shows up. This happens, for 
example, at the step between Fi and F2 for Fl = Fr, as 
seen in Fig.|31 The peak height, 0.5625, exceeds that of 
the adjacent plateaus, 0.5 and 5/9. These features can 
even appear in the regime of negligible Coulomb charging 
energy, as previously shown in Ref. 'I8'. However, we 
point out that the behavior at the steps is going to be 
strongly affected by second-order tunneling events, so- 
called cotunneling. For example, the current steps will 
show additional broadening due to the intrinsic linewidth 
F. On the other hand, the plateau values in regimes 1 - 
4 will not be affected by second-order effects. 

In summary, we presented a theory of transport 
through a molecule or quantum dot that allows for a sys- 
tematic perturbation expansion of current and shot noise 
in the coupling strength between molecular orbital (dot 
level) and electrodes. We analyzed first-order transport 
in detail. In a complete overview, we derived analytic 
expressions for current, noise, and the Fano factor for all 
possible situations involving a single interacting level. In 
particular, we discussed how the sequence of Fano factor 
plateaus can provide information about the asymmetry 
ratio of the coupling strengths as well as the presence of 
Coulomb charging energy. 
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